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Abstract 



We show that the initial-state interactions from gluon exchange between the 
incoming quark and the target spectator system lead to leading-twist single-spin 
asymmetries in the Drell-Yan process -ffi-ffj ~^ £ + £~X. The QCD initial-state 
interactions produce a T— odd spin-correlation Sh 2 'Ph x x Q between the target 
spin and the virtual photon production plane which is not power-law suppressed 
in the Drell-Yan scaling limit at large photon virtuality Q 2 at fixed xf- The 
single-spin asymmetry which arises from the initial-state interactions is not re- 
lated to the target or projectile transversity distribution 5qH(x, Q). The origin 
of the single-spin asymmetry in np^ — > l + £~X is a phase difference between two 
amplitudes coupling the proton target with J? = ±| to the same final-state, the 
same amplitudes which are necessary to produce a nonzero proton anomalous 
magnetic moment. The calculation requires the overlap of target light-front 
wavefunctions with different orbital angular momentum: AL Z = 1; thus the 
SSA in the Drell-Yan reaction provides a direct measure of orbital angular mo- 
mentum in the QCD bound state. The single-spin asymmetry predicted for the 
Drell-Yan process trp^ — ► £ + £~X is similar to the single-spin asymmetries in 
deep inelastic semi-inclusive leptoproduction £p' — ► £'ttX which arises from the 
final-state rescattering of the outgoing quark. The Bjorken-scaling single-spin 
asymmetries predicted for the Drell-Yan and leptoproduction processes high- 
light the importance of initial- and final-state interactions for QCD observables. 



1 Introduction 

Single-spin asymmetries in hadronic reactions provide a remarkable window to QCD 
mechanisms at the amplitude level. In general, single-spin asymmetries measure 
the correlation of the spin projection of a hadron with a production or scattering 
plane fl[]. Such correlations are odd under time reversal, and thus they can arise in a 
time-reversal invariant theory only when there is a phase difference between different 
spin amplitudes. Specifically, a nonzero correlation of the proton spin normal to a 
production plane measures the phase difference between two amplitudes coupling the 
proton target with J* = ±| to the same final-state. The calculation requires the 
overlap of target light-front wavefunctions with different orbital angular momentum: 
AL Z = 1; thus a single-spin asymmetry (SSA) provides a direct measure of orbital 
angular momentum in the QCD bound state. 

Consider the SSA produced in semi- inclusive deep inelastic scattering tp^ — > i'nX. 
In the target rest frame, such a single target spin correlation corresponds to the T-odd 
triple product S p ■ p n x q. (The covariant form of this correlation is e tlvaT S£'p 1 ' 'q° p^-) 
Significant asymmetries Ajjl and Ay? of this type have in fact been observed for 
targets polarized parallel to or transverse to the lepton beam direction 0, |3|. 

In a recent paper 0] we have shown that the QCD final-state interactions (gluon 
exchange) between the struck quark and the proton spectator system in semi-inclusive 
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deep inelastic lepton scattering can produce single-spin asymmetries which survive in 
the Bjorken limit. Such effects are proportional to the matrix element of a higher- 
twist quark-quark-gluon correlator in the target hadron, and thus it has been assumed 
on dimensional grounds that any SSA arising from this source must be suppressed 
by a power of the momentum transfer Q in the Bjorken limit. However, another 
momentum scale enters into the semi-inclusive process — the transverse momentum 
t± = — q± of the emitted pion relative to the photon direction, and we have shown 
that the power-law suppression due to the higher-twist quark-quark-gluon correlator 
takes the form of an inverse power of r± rather Q. As shown in the Appendix, r\ can 
be written in terms of the invariant momentum transfer squared t from the proton to 
the spectator system and the Bjorken variable. 

Corrections from spin-one gluon exchange in the initial- or final-state of QCD 
processes are not suppressed at high energies because the coupling is vector-like. 
Therefore, as a consequence of the gauge coupling of QCD, single-spin asymmetries 
in semi-inclusive deep inelastic scattering survive in the Bjorken limit of large Q 2 
at fixed Xbj and fixed r_j_. Recently it has been shown || [| that the same type 
of final-state interaction is the origin of the leading-twist diflractive component in 
deep inelastic scattering, implying that the pomeron is not a universal property of 
the target proton's wavefunction, and that it depends in detail on the deep inelastic 
scattering (DIS) process itself. Diflractive processes in DIS in turn lead to nuclear 
shadowing in the case of nuclear targets, showing that shadowing is not an intrinsic 
property of nuclear wavef unctions. 

The final-state phases which we compute are analogous to the "Coulomb" phases 
to the hard subprocess which arises from gauge interactions between outgoing charge 
particles in QED 0. More specifically, we require the difference between the gauge 
interaction phases for the J* = ±| amplitudes. The phases depend on the spin 
because the outgoing particles interact at different impact separation corresponding 
to their different relative orbital angular momentum. 

In our previous paper EL we explicitly evaluated the SSA for electroproduction 
for a specific model of a spin-| proton of mass M with charged spin-| and spin-0 
constituents of mass m and A, respectively, as in the QCD-motivated quark-diquark 
model of a nucleon. The basic leptoproduction reaction is then 7*p — > q(qq)o- Our 
analysis predicts a nonzero SSA for the target spin normal to the photon to quark- 
jet S p ■ p q x q which can be determined by using a jet variable such as thrust to 
determine the current quark direction; i.e., we predict a SSA even without final-state 
jet hadronization. Our mechanism is thus distinct from a description of SSA based 
on transversity and phased fragmentation functions. 

Recently Collins |§ has pointed out some important consequences of these results 
for SSA in deep inelastic scattering. In his treatment the final-state interactions of 
the struck quark are incorporated into Wilson line path-ordered exponentials which 
augment the light-cone wavefunctions. [See also 0.] Since the final-state inter- 
actions appear at short light-cone times Ax + = 0{l/v) after the virtual photon 
acts, they can be distinguished from hadronization processes which occur over long 
times. Collins has stressed the fact that single-spin asymmetries probe the partonic 
structure associated with chiral-symmetry breaking. Furthermore, these results show 
that time-reversal-odd parton densities are allowed, opening up a whole range of 
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Figure 1: The initial-state interaction in the Drell-Yan process. 



phenomenological applications |T0| , [TT , [12| . In particular, as noted by Collins, initial- 
state interactions between the annihilating antiquark and the spectator system of the 
target can produce single-spin asymmetries in the Drell-Yan process. 

In this paper we shall extend our analysis to initial- and final-state QCD effects to 
predict single-spin asymmetries in hadron-induced hard QCD processes. Specifically, 
we shall consider the Drell-Yan (DY) type reactions |1| such as ixp^ — > £ + £~X. Here 
the target particle is polarized normal to the pion production plane. The target 
spin asymmetry can be produced due to the initial-state gluon-exchange interactions 
between the interacting antiquark coming from one hadronic system and the spectator 
system of the other. This is shown in the diagram of Fig. 1. The importance of initial- 
state interactions in the theory of massive lepton pair production, Q± broadening, 
and energy loss in a nuclear target has been discussed in Refs. [[b|, [L5|, |T6|j . 

The orientation of the target spin S z = ±1/2 corresponds to amplitudes differing 
by relative orbital angular momentum AL Z = 1. The initial-state interaction from 
a gluon exchanged between the annihilating antiquark and target spectator system 
depends in detail on this relative orbital angular momentum. In contrast, the initial 
or final-state interactions due to the exchange of gauge particles between partons not 
participating in the hard subprocess do not contribute to the SSA. Such spectator- 
spectator interactions occur at large impact separation and are not sensitive to just 
one unit difference AL Z = 1 of the orbital angular momentum of the target wave- 
function. 

Our mechanism thus depends on the interference of different amplitudes arising 
from the target hadron's wavefunction and is distinct from probabilistic measures 
of the target such as transversity. It is also important to note that the target spin 
asymmetries which we compute in the DY and DIS processes require the same over- 
lap of wavefunctions which enters the computation of the target nucleon's magnetic 
moment. In addition, by selecting different initial mesons in the DY process, we can 
isolate the flavor of the annihilating quark and antiquark. The flavor dependence of 
single-spin asymmetries thus has the potential to provide detailed information of the 
spin and flavor content of nucleons at the amplitude level. 

As in our analysis of semi- inclusive DIS, we shall calculate the single-spin asym- 
metry in the Drell-Yan process induced by initial-state interactions by adopting an 
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effective theory of a spin-^ proton of mass M with charged spin-| and spin-0 con- 
stituents of mass m and A, respectively, as in a quark-diquark model. We will take 
the initial particle to be just an antiquark. The result for specific meson projectiles 
such as Mp^ — > £ + £~X is then obtained by convolution with the antiquark distribu- 
tion of the incoming meson. One can also incorporate target nucleon wavefunctions 
with a quark-vector diquark structure. In a more complete study, one should allow 
for a many-parton light-front Fock state wavefunction representation of the target. 
The results, however, are always normalized to the quark contribution to the proton 
anomalous moment, and thus are basically model- independent. 



2 Crossing 

There is a simple diagrammatic connection between the amplitude describing the 
initial-state interaction of the annihilating antiquark, which gives a single-spin asym- 
metry for the Drell-Yan process Tip 1 * — > £ + £~X, and the final-state rescattering am- 
plitude of the struck quark, which gives the single-spin asymmetries in semi-inclusive 
deep inelastic leptoproduction £p^ — > fnX. The crossing of the Feynman amplitude 
for 7*(g)p(-P) — > {q + r)(P — r) in DIS gives (— q — r)p(P) — > 7*(— q)(P — r) for DY 
by reversing the four- vectors of the photon and quark lines. The outgoing quark with 
momentum q + r in DIS becomes the incoming antiquark with momentum —q — r 
in DY. We can use crossing of the Lorentz invariant amplitudes for DIS as a guide 
for obtaining the amplitudes for DY amplitude |]17|| . [In the next section it will be 
convenient to label q = —q with q + > 0. ] 

In general, one cannot use crossing to relate imaginary parts of amplitudes to each 
other, since under crossing, real and imaginary parts become connected. However, in 
our case, the relevant one-gluon exchange diagrams in DIS and DY are both purely 
imaginary at high energy, so their magnitudes are related by crossing. Thus a crucial 
test of our mechanism is an exact relation between the magnitude and flavor depen- 
dence of the SSA in the Drell-Yan reaction and the SSA in deep inelastic scattering. 
We thus predict the DY SSA of the proton spin with the normal to the antiquark to 
virtual photon plane: S p -pgX q. It is identical - up to a sign - to the SSA computed 
in DIS for S p ■ p q x q. 

The phase arising from the initial- and final-state interactions in QCD is analogous 
to the Coulomb phase of Abelian QED amplitudes. The Coulomb phase depends on 
the product of charges and relative velocity of each ingoing and outgoing charged 
pair Q. Thus the sign of the phase in DY and DIS are opposite because of the 
different color charge of the ingoing 3c antiquark in DY and the outgoing 3c quark 
in DIS. In order to check the sign, we will carry out the DY calculation explicitly in 
the next section. 

The asymmetry in the Drell-Yan process is thus the same as that obtained in DIS, 
with the appropriate identification of variables, but with the opposite sign. This has 
been stressed recently by Collins || . Therefore the single-spin asymmetry transverse 
to the production plane in the Drell-Yan process can be obtained from the results of 
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our recent paper 
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Here A = = where z/ is the energy of the lepton pair in the target rest frame. 
The kinematics are given in detail in the Appendix. 



3 Calculation 



We can check the results of the previous section obtained using crossing by performing 
a direct calculation of the qp^ — > 7*(gg)o amplitude where we take a spin- zero diquark 
for the proton spectator system. The kinematics are (q — r)p(P) — > r j*{q){P — r), as 
in Fig. 1. The J z = +| two-particle Fock state is given by []I8J |T^] 



two particle V ' -L 

d 2 ki dx 



where 
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(2) 



^ +1 _(x,k ± ) = (M+^)<p , 



(3) 



The scalar part of the wavefunction ip depends on the dynamics. In the perturbative 
theory it is simply 



if = (f{x, k ± ) 



'l-x 



M 2 



l-x 



(4) 



In general one normalizes the Fock state to unit probability. 

Similarly, the J z = — | two-particle Fock state has two components 

ij l _l(x,k ± ) = (M+^)<p . 



(5) 



The spin-flip amplitudes in (|3|) and (|5|) have orbital angular momentum projection 
l z = +1 and —1 respectively. The numerator structure of the wavefunctions is char- 
acteristic of the orbital angular momentum, and holds for both perturbative and 
non-perturbative couplings. 



6 



We require the interference between the tree amplitude of Fig. 1(a) and the one- 
loop amplitude of Fig. 1(b). The contributing amplitudes have the following structure 
through one-loop order: 

: ( M +^)C(h + i e -^ gi ) (6) 

- (±^) C<* + ^) (7) 

= (^W + i ^ S2 ) (8) 

= (M + ^)C(h + i^ gi ), (9) 

where 





H) 




H) 




♦T) 







C = - 5ei P + ^-A2A(1-A) (10) 

1 

ri + A(l- A)(~M 2 + ^ + ^ 



^ = ^ , a„ ^ ^ , , A 2 i • ( U ) 



The label -fT / 4J- corresponds to J* = ±~ of the proton spin. The second label | / I 
gives the spin projection J| = ±| of the interacting spin-| constituent antiquark of 
the other proton. Here e% and e 2 are the electric charges of the proton's constituents 
q and (qq)o, respectively, and g is the coupling constant of the effective proton-g-(gg) 
vertex. The first term in (Q) to (|9|) is the Born contribution of the tree graph. The 
crucial result will be the fact that the contributions g\ and g 2 from the one-loop 
diagram Fig. 1(b) are different, and that their difference is infrared finite. A gauge 
boson mass X g will be used as an infrared regulator in the calculation of g\ and g 2 • The 
final result for g\ — g 2 is infrared finite, and X g can be set to zero. The calculation 
will be done using light-cone time-ordered perturbation theory, or equivalently, by 
integrating Feynman loop-diagrams over dk~ . 

We take q to lie in the z — x plane, q = (q x , q y , q z ) = (q 1 , 0, q 3 ). We denote q + as 

q + = j3 P+ . (12) 



From energy conservation, we get 



9 ~ - pmFa] ' (13) 



Then we have the relation 



ft = ft<T ~ ft = ^ ft , (14) 

where 

P > A . (15) 
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Further details on the kinematics are given in the Appendix. 

The covariant expression for the four one-loop amplitudes of diagram Fig. 1(b) 

is: 



A oae - loo »(l) (16) 
d A k 



ig e 1 { -e x e 2 
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d 2 k ± r Af(I) 



-i 9 e 1 (-e 1 e 2 ) J ^ J P + dx 



x dk 



2(2tt) 4 7 P+ 4 x (j3 — x) (x - A) (1 - x) 

1 
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1 

where we used A; 4 " = xP + . The numerators A/*(7) are given by 

TH 

mt^i) = N (M + -) (17) 

x 

m^i) = n { +kl - ik2 ) (is) 

X 

— k 1 — ilc 2 

m^) = n ( — ) (19) 

777 

AA(^t) = ^ (M H ) , (20) 

x 

where 



AT = 2P + Jf3-A x q- (-P + [(l -x) + (1 - A)]) , (21) 



and q~ = P +^_ A ^ as given in (J13|) . For the [current]- [gauge boson propagator]- 

[current] factor, in Feynman gauge only the —g + ~ term of the gauge boson propagator 
— g^ v contributes in the Bjorken limit, and it provides a factor proportional to q~ in 
the numerator which cancels the q~ in the denominator which provides the imaginary 
part. Therefore the result scales in the Bjorken limit. 

The integration over k~ in (O) does not give zero only if < x < 1. We first 
consider the region A < x < 1. 

^one-loop (j) (22) 

d 2 k ± r . , A/"(J) 



/tfc ft 1 /* 
2 / 27r )4 y p+rfx 



2(27r) 4 7 P+ 4 x (ft — x) (x — A) (1 - x) 
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r (i-x)p+ xp+ j I (i-x)p+ y ip-x)p+ 
1 



_ (A2+fcj_)-, e _ _ (Ag+(fc x -r ± )2)-ie 
^ (l-x)P+ ' (x-A)P+ 



where we used k + = xP + . The result is identical to that obtained from light-cone 
time-ordered perturbation theory. 

The phases Xi needed for single-spin asymmetries come from the imaginary part 
of (0), which arises from the potentially real intermediate state allowed before the 
rescattering. The imaginary part of the propagator (light-cone energy denominator) 
gives 

- , , P - M (A 2 + £i) , _ (m 2 + (-k ± + q ± ) 
— mo —P + 



[l-x)P+ * (J3-x)P+ 

J_ (g - A) 2 
P+ q\ 

where 



™ =5 — A — d>) , (23) 



? = 2(^-A) gL ' ( ^" f - L) . (24) 

<?_i_ 

Since the exchanged momentum 5P + is small, the light-cone energy denominator 

corresponding to the gauge boson propagator is dominated by the ^ k± ^r^r— ~ term. 

This gets multiplied by (x — A), so only (k± — r±) 2 + A 2 , appears in the propagator, 
independent of whether the photon is absorbed or emitted. The contribution from 
the region < x < A thus compliments the contribution from the region A < x < 1. 

We can integrate ([22]) over the transverse momentum using a Feynman parame- 
trization to obtain the one-loop terms in (M) to (pi). 



»i 

■I ' ' /_ rfo — — : x2 , . , 7T7, , , , > . T— (2- r )) 



o ' " a (l - «)rl + aA 2 , + (1 - a )A(l - A)(-M 2 + + ^ 



92 Jo ^ a{l - a)rl + a\ 2 g + (1 - a)A(l - A)(-M 2 + ^ + ^) ' (26) 
We define: 

V z = C- 1 ( |A(^T)| 2 - W^T)| 2 + \A(f^i)\ 2 - \A(^l)\ 2 ) (27) 

= c- 1 ( ( A(ft-T)M(^-T) + A(ft-T)^-T)* ) (28) 
+ ( A(t^iyA(^l) + ) ) 

V y = (T 1 ( i ( A(1Ht)M(^t) - AflHTMGHT)* ) (29) 



where the normalization from the unpolarized cross section is 



C = \A(^)\ 2 + |A(^T)| 2 + \A(^i)\ 2 + |A(^I)| 2 . (30) 

We can assume for convenience that the initial-state interactions generate a phase 
when exponentiated, as in the Coulomb phase analysis of QED. The rescattering 
phases e %Xi (i = 1,2) with \% — tan _1 (^^f-) are thus distinct for the spin-parallel 
and spin-antiparallel amplitudes. The difference in phase arises from the orbital 
angular momentum k± factor in the spin-flip amplitude, which after integration gives 
the extra factor of the Feynman parameter a in the numerator of g% compared to g\, as 
we can see in ( ^5|) and ( f26|) . Notice that the phases Xi are each infrared divergent for 
zero gauge boson mass X g — > 0, as is characteristic of Coulomb phases. However, the 
difference Xi ~ X2 which contributes to the single-spin asymmetry is infrared finite. 
We have verified that the Feynman gauge result is also obtained in the light-cone 
gauge using the principal value prescription. The small numerator coupling of the 
light-cone gauge boson is compensated by the small value for the exchanged l + = 5P + 
momentum. 

The virtual photon and produced hadron define the production plane which we 
will take as the z — x plane. From Eqs. and (|29|), the azimuthal single-spin 

asymmetry transverse to the production plane is given by 
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in agreement with the crossing properties described in section 2. The linear factor of 
r 1 = r x reflects the fact that the single-spin asymmetry is proportional to S p -qxr since 

q ~ —z\q\ and S p = rky. The kinematics are given in more detail in the Appendix. 
The prediction for V y as a function of A and r± is identical but with opposite sign 
to that illustrated in Fig. 4 of Ref. [|]. 

Our analysis can be generalized to the corresponding calculation in QCD. The 
initial-state interaction from gluon exchange has the strength — > CVa s (/z 2 ). The 
scale of a s in the MS scheme can be identified with the momentum transfer carried 
by the gluon /i 2 = e~ 5//3 (fcj_ — f±) 2 |p0 |. The matrix elements coupling the proton 



to its constituents will have the same numerator structure as the perturbative model 
since they are determined by orbital angular momentum constraints. The strengths 
of the proton matrix elements can be normalized quark by quark according to their 
contributions to the target nucleon's anomalous magnetic moment weighted by the 
quark charge squared. In QCD, r± is the magnitude of the momentum of the cur- 
rent quark jet relative to the virtual photon direction. Notice that for large r±, V y 

decreases as aa ^ r A) Xb ^ Iri - lnr ± _ physical proton mass M appears since it is present 

in the ratio of the L z = 1 and L z = matrix elements. This form is expected to be 
essentially universal. 
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4 Summary 



We have shown that the same physical mechanism which produces a leading-twist 
single-spin asymmetry in semi- inclusive DIS, also leads to a leading twist single-spin 
asymmetry in the Drell-Yan process. The initial-state interaction between the an- 
nihilating antiquark with the spectator of the target produces the required phase 
correlation. The equality in magnitude, but opposite sign, of the single-spin asymme- 
tries in semi-inclusive DIS and the corresponding Drell-Yan processes is an important 
check of our mechanism. 

It has been conventional to assume that the effects of initial- and final-state in- 
teractions are always power-law suppressed for hard processes in QCD. In fact, this 
is not in general correct, as can be seen from our analyses of leading-twist single-spin 
asymmetries in the Drell-Yan process and semi-inclusive deep inelastic scattering. 
The initial- and final-state interactions which survive in the scaling limit occur in 
light-cone time r = 0(1/ Q) immediately before or after the hard subprocess. Other 
initial- and final-state interactions, such as those between the spectator of the incident 
hadron and the spectator of the target hadron in the DY process, take place over long 
time scales, and they only provide inconsequential unitary phase corrections to the 
process. This is in accord with our intuition that interactions which occur at distant 
times cannot affect the primary reaction. 

A natural framework for the wavefunctions which appear in the SSA calculations 



is the light-front Fock expansion [23]. In principle, the light-front wavefunctions 
for hadrons can be obtained by solving for the eigen-solutions of the light-front QCD 
Hamiltonian. Such wavefunctions are real and include all interactions up to a given 
light-front time. The final-state gluon-exchange corrections which provide the SSA 
for semi-inclusive DIS occurs immediately after the virtual photon strikes the ac- 
tive quark. Such interactions are not included in the light-front wavefunctions, just 
as Coulomb final-state interactions are not included in the Schrodinger bound state 
wavefunctions in QED. Collins || has argued that since the relevant rescattering 
interactions of the struck quark occur very close in light-cone time to the hard inter- 
action, one can augment the light-front wavefunctions by a Wilson line factor which 
incorporates the effects of the final-state interactions in semi- inclusive DIS. However, 
such augmented wavefunctions are not universal and process independent; for exam- 
ple, in the case of the DY process, an incoming Wilson line of opposite phase must 
be used. 

Our formalism can be adopted to single-spin asymmetries in more general hard 
inclusive reactions, such as pp^ — > 7rX, where the pion is detected at high transverse 
momentum |23|, p4 |. In such cases one must identify the hard quark-gluon subprocess 



and analyze a set of gluon exchange corrections which connect the spectators of the 
polarized hadron with the active quarks and gluons of the hard subprocess. An ex- 
ample of a final-state interaction which can cause a single asymmetry in pp" — > ttX is 
shown in Fig. 2. However, this type of final-state interaction cannot be readily iden- 
tified as an augmented target wavefunction. It is also clear from our analyses that 
there are potentially important corrections to the hard quark propagator in hard 
exclusive subprocesses such as deeply virtual Compton scattering or exclusive me- 
son electroproduction. These rescattering interactions of the propagating quark can 
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Figure 2: An example of a final-state interaction which can cause a single asymmetry 
in pp^ — ► ix X . 



provide new single-spin observables and will correct analyses based on the handbag 
approximation. 

It should be emphasized that the same overlap of light-front wavefunctions with 
AL Z = 1 which gives single-spin asymmetries also yields the Pauli form factor F 2 (t) 
and the generalized parton distribution E(x,£,t) entering deeply virtual Compton 



scattering [p5| , |26| , |27| , Pq , P9| . Each quark of the target wavefunction appears addi- 
tively, weighted linearly by the quark charge in the case of the Pauli form factor and 
weighted quadratically in the case of deep inelastic scattering, the Drell-Yan reaction 
and deeply virtual Compton scattering. 

The empirical study of single-spin asymmetries in hard inclusive and exclusive 
processes thus provides a new window to the investigation of hadron spin, angular 
momentum, and flavor structure of hadrons. 
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Appendix: Kinematics 

From Fig. 1 we have [a = (a + , a~, a±)] 

P = (p+ , ^ , ± ) (32) 

p ' = ( (1 - A)F+ - (T^- - rl ) (33) 

P~, = (VI - A)P + , m \+ { _%;} )2 , U - fx) , (35) 
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where 

s - m 2 - M 2 = 2P r P = m2+ ^2) r±)2 + M ' (/? " A) ' (36) 



Energy conservation gives 

m 2 + (q ± -f ± ) 2 2 _ g 2 + qj A 2 + rj 

H3-A) +M ~ f3 + (1-A)- (37) 

Here s and q 2 are large. It is convenient to work in a frame where (3 — A = 0(1) and 
q\ is large so that s ~ wza- In such a frame, q ~ — z\q\. Then Eq. d37|) gives 



hi 0-A 



From (H), (|3|) and (g§ we have 



2P.g = ^ + /5M 2 (39) 



and thus to leading twist 



IP -q ' 

where P ■ q/M = u is the energy of the lepton pair in the target rest frame. 
From (|32|) and (B3|) we have 



(40) 



2P-P' = + (1 - A)M 2 . (41) 

Since t = (P - P') 2 , we have 2P • P' = -t + M 2 + A 2 . Therefore, (0) gives 

= (1 - A)(-t + M 2 + A 2 - (1 - A)M 2 ) - A 2 . (42) 
which relates to invariants. 
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